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1. Introduction 

The advent of direct numerical simulations of turbu- 
lence has opened new avenues for research on turbu- 
lence physics and turbulence modeling. Although lim- 
ited to relatively low Reynolds numbers, direct sim- 
ulations have already provided new insight into the 
structure of turbulence. Direct numerical simulation 
provides values for anything that the scientist or mod- 
eler would like to know about the flow, and this is be- 
ginning to have tremendous payoff. Simulations have 
enabled detailed evaluations to be made of turbulence 
models intended for use in engineering codes, as well as 
models of the small-scale turbulence required for Large 
Eddy Simulations (LES), a promising technique for 
practical numerical simulations of turbulence at high 
Reynolds numbers. This paper presents an overview 
of some recent advances in the physical understand- 
ing of turbulence and in turbulence modeling obtained 
through such simulations. 

One class of turbulence that has been studied exten- 
sively by direct simulations is homogeneous turbulence . 
In homogeneous turbulence, the mean velocity gradi- 
ents are prescribed and must be uniform in space, and 
all turbulence statistics are independent of position. 
One of the difficult problems in turbulence simulations 
is the provision of boundary conditions. It has been 
demonstrated that, in simulating homogeneous flows, 
one may use periodic boundary conditions, provided 
that the computational domain is larger than twice the 
distance between points where statistical motions are 
significantly correlated. Hence, the boundary motions 
do not have to be known to simulate homogeneous tur- 
bulence, but instead emerge as a result of the dynamics 
of the flow and the periodicity. It is for this reason that 
much of the first work on direct simulations has dealt 
with homogeneous turbulence. 

A program developed by Dr. Robert Rogallo (1981) 
at Ames has become the bulwark of direct simula- 
tions of various homogeneous turbulent flows. The 
program uses de-aliased spectral methods to obtain 
high accuracy at small scales. It works in a coordinate 
system that deforms with the mean flow using vari- 
ables that are computed exactly in the limit of very 
rapid distortions, and so the program has been used 
to cover a wide range of mean deformation conditions. 
The Rogallo code has been used to study the decay of 
isotropic turbulence, the response of isotropic turbu- 
lence to imposed mean strain and rotation, and homo- 


geneous shear flow, including scalar transport. This 
work spans the Ames supercomputers from ILLIAC- 
IV to CRAY-2, and has recently been adapted for use 
on Masscomp supermini computers. 

This paper highlights the results of three recent Stan- 
ford Ph.D. Dissertations, carried out on the Cray- 
XMP at NASA/Ames under the direction of the au- 
thor as part of a cooperative program of turbulence re- 
search with colleagues at Ames. This program, which 
has spanned well over a decade, has produced many of 
the currently active researchers in this field, and has 
been instrumental in the establishment of Ames as the 
world center for this type of research. A key outgrowth 
of this prorgam is the recently-formed NASA/Stanford 
Center for Turbulence Research , which will bring to- 
gether experimentalists, computational scientists, and 
modelers in a concerted attack on the physics and 
modeling of turbulence. 

Two of the disserations discussed below deal with ho- 
mogeneous turbulence, and were carried out using the 
Rogallo code. They represent the current state-of-the 
art in simulations of such flows, and illustrate the sort 
of insight into turbulence physics and models that can 
be obtained from direct simulations. The third disser- 
tation was the first direct simulation of the spatially- 
developing turbulent mixing layer, and is characteris- 
tic of the state of affairs in the direct simulation of 
flows of practical interest. Current work in extending 
this simulation to the supersonic mixing layer, for ap- 
plications in both aerodynamics and scramjet engines, 
will be mentioned. 

2. Response of turbulence to irrotational strain 
Two of the building-block flows for turbulence mod- 
elers are (l) homogeneous turbulence undergoing in- 
compressible, irrotational mean strain and (2) relax- 
ation to isotropy after cessation of the mean strain- 
ing. These two situations apply approximately to the 
turbulence (l) in the contraction section of low-speed 
wind tunnels and (2) in the straight test section fol- 
lowing the nozzle. Hence, this idealized case has prac- 
tical importance, especially in wind tunnels designed 
for low free-stream turbulence. 

Different types of mean deformation are of inter- 
est. Wind tunnels are either axisymmetric or two- 
dimensional, and the behavior of the turbulence is 
quite different in these cases. A simple model based on 
the idea that the turbulence consists of a complex tan- 
gle of vortex filaments is useful in understanding these 
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differences. The streamwise strain stretches the vortex 
filaments that form the turbulence, concentrating the 
vorticity and increasing the intensity of the swirling 
motion around each filament. The cross-stream con- 
traction pushes these filaments closer together. In ax- 
isymmetric contraction the filaments are pushed to- 
gether from both sides, whereas in plane strain (e.g. 
two-dimensional mean contraction), the vortex fila- 
ments are pushed together only in one direction. 

In the diffuser segment of a wind tunnel, a still dif- 
ferent type of mean deformation is imposed. In a 
round diffuser, vortex filaments are stretched in the 
two cross-stream directions, and compressed in the 
streamwise direction. This deforms the filaments into 
thin sheets, giving rise to intense local shear layers in 
the flow. 

This line of thinking suggests that the small-scale 
structure of the turbulence, as determined by the tan- 
gle of vortex filaments, should be different under differ- 
ent mean distortions. But one of the common premises 
in turbulence modeling is that the small-scale struc- 
ture is isotropic, which is not consistent with the model 
of distorted filaments. The model also suggests that 
the rate of return to isotropy upon cessation of mean 
deformation should depend on the nature of the vortex 
filament tangle, i.e. on the nature of the prior mean 
strain. These are features that one would expect to 
see in a numerical simulation and would like to see 
exhibited by analytical model of turbulence. 

Moon Lee (1986), now a Postdoctoral Researcher at 
NASA/Ames, performed an extensive series of direct 
simulations of homogeneous turbulence subjected to 
a variety of irrotational strains, and also studied the 
relaxation upon cessation of mean deformation. The 
objective of his dissertation was to use the simulation 
data to evaluate second-order closure models of tur- 
bulence. His simulations used a 128 3 mesh, a second- 
order Runga-Kutta time advance, carefully-configured 
initial conditions, and were executed at Reynolds num- 
bers, based on the length and velocity scales of the 
energy-containing motions, of the order of 50. Some 

background necessary for understanding the impact of 
his work will now be outlined. 

In second-order closures, one deals with the Reynolds 
stress tensor 

= W (!) 

or, alternatively, with the turbulent velocity scale 

9 2 = Rii (2a) 

and the Reynolds stress anisotropy tensor 

% = (Rii - p 2 S i} )/q\ (26) 

In isotropic turbulence 6*y = 0. If the energy in one 
component, say u^, goes to zero, then 6n = — and 


if all of the fluctuation energy is in one component, say 
u^, then &n = |. Thus, the values of 6*y are limited. 
These limits are conveniently described in terms of the 
invariants of the 6*y tensor. The first invariant is /*> = 
bu = 0. The second and third invariants are 

lib = ~~ ^bijbji (3a) 

and 

II h = - bijbjkbki . (36) 

The limits outlined above translate into limits on Ih 
and Illbj and one finds that all possible states of tur- 
bulence must lie within a certain area on the Illb - 
lib plane, which is called the anisotropy invariant map 
(AIM) for the Reynolds stress. Simple eddy viscosity 
models for the Reynolds stresses do not necesarily sat- 
isfy this condition, particularly in cases of very strong 
strain-rate, and hence sometimes these models pro- 
duce unrealizable turbulence. The AIM is therefore a 
very useful tool for assessing turbulence models. 

Second-order closure models for homogeneous flows 
require closures for the equations of evolution of i^y 
and e , the rate of dissipation of mechanical energy 
per unit volume by the turbulence. Alternatively, the 
evolution equations for 6 t y and q 2 are used in place of 
those for R,y . The turbulent kinetic energy equation 
describes the evolution of q 2 } and is 

% = (4«) 

where 

P = -RijSi y (46) 

is the rate of energy transfer from the mean motion to 
the turbulence ( production ) and 

Sii = \{U^ + U s%i ) (4c) 

is the mean strain-rate (commas denote partial differ- 
entiation). No modeling is necessary in (4a) because it 
contains only the model variables and the prescribed 
mean velocity gradients. 

The equation for 6*y is 

= Aij — ^2 (^»y — 26*y). (5) 

Here A t y is a production term involving only the an- 
isotropy tensor b and the mean velocity gradients, 
and hence does not require modeling. <£ t y, which in- 
volves the anisotropy of the viscous dissipation and the 
pressure-strain terms, must be modeled. It contains a 
rapid pressure- strain part that changes instantly when 
the mean velocity gradients are changed, which must 
be modeled as being proportional to the mean gra- 
dients, and slow terms that do not contain the mean 
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velocity gradients. For the special case of no mean gra- 
dients (the return to isotropy problem), contains 
only the slow terms. Hence, the modeling of these 
terms can be examined by studying direct simulations 
of the return to isotropy, and this was the primary 
objective of Moon Lee’s dissertation. 

If one is going to work with bij , q 2 > and e as variables, 
then the model for <f>ij must be expressed interms of 
these variables. The tensor character must be car- 
ried by tensors developed from bij . By the Cayley- 
Hamilton theorem, there are only two linearly inde- 
pendent tensors that can be involved. Thus, the most 
general model of the slow part of </>ij possible in this 
type of second-order closure is 


The vorticity anisotropy tensor is defined by 


V*i-fr a 6ij 

w 2 


(8d) 


and its invariants, //„ and II I v are defined in the same 
way as those for b. 

The dissipation is 


Dij = 2 


(9a) 


The trace of this tensor is 

Du = 2e (96) 


4>ij = (<* + 2)bij + ft(b ik b k j H- -II b 8ij) (6) 

where a and ft are scalar coefficients that can depend 
upon the scalar invariants II b and III b and on other 
scalars, such as the turbulent Reynolds number Rp = 
q A /(ei/). An immediate consequence of this model is 
that the return to isotropy must follow a trajectory on 
the Reynolds stress AIM given by 

dllh = 3a/// b + |/?// b 2 

dll b 2 ocII b - ZpIII b ‘ [ 1 

Hence, the trajectory should be a unique function of 
the position on the AIM (and i?^). Thus, by following 
the trajectory in the simulations, one should be able 
to deduce the parameters a and ft, provided that the 
assumed model is adequate to describe the flow. 

The simulations showed that the return trajectory 
is no t a unique function of the AIM state; the most 
striking deviations were in the case of axisymmetric ex- 
pansion flows, where turbulence that had been rapidly 
strained actually moved further away from isotropy 
when the strain-rates was removed, while turbulence 
that had been strained slowly to the same anisotropy 
moved towards isotropy upon the removal. Thus, at 
least for the case of strong strain-rate (but of the or- 
der of the strain rate in many practical situations), the 
standard second-order model simply is inadequate. 

We then turned to looking for alternative ways to 
characterize and model the turbulence, using two other 
tensors characterizing the anisotropy of the small-scale 
motions. Analogous to q 2 and 6 t *y we define the 
vorticity tensor by 

V i} - = w'w'- (8a) 

where wj is a turbulent vorticity component. The 
mean-square vorticity is 

w 2 = Vu (86) 

which, for homogeneous turbulence, is related to the 
dissipation by 

e = i/a; 2 . (8c) 


The associated dissipation anisotropy tensor is defined 
by 


dii = 


Dij fefty 
2e 


(9c) 


and its invariants, lid and II Id are defined in the same 
way as those for b. 

The AIMs for v and d look exactly like that for b, 
although the phisical meanings of the limiting bound- 
aries are different. For example, turbulence in which 
the vortex filaments have been stretched infinitely in 
one direction has a one-dimensional vorticity field but 
a two-dimensional velocity field, and hence its state is 
at different points on the v and b AIMs. 

The simulations revealed a very surprising result that 
shakes the roots of turbulence theory. The anisotropy 
of the small-scale turbulence was not found to be much 
smaller than the anisotropy of the large-scale turbu- 
lence, as the lore of turbulence suggests. Instead, 
there tended to be strong relatonships between the 
anisotropies at small and larges scales, as revealed by 
the instantaneous states of the turbulence on the vor- 
ticity /dissipation AIMS and the Reynolds stress AIM, 
respectively. It has been suggested that this is a re- 
sult of the low Reynolds numbers used in the simu- 
lations. However, we can see no trend for the small- 
scale anisotropy to be less at higher Reynolds numbers. 
Figure 1 shows the second invariants of the vortic- 
ity and Reynolds stress anisotropies during the relax- 
ation from simple strains. Note that for virtually ev- 
ery field examined the vorticity anisotropy was larger 
than the Reynolds stress anisotropy, as measured by 
their second invariants. Figure 2 shows the dissipation 
anisotropy vs. that for the Reynolds stress for these 
same velocity fields. Note that the second invariants of 
these anisotropies are generally of comparable value, 
again indicating the non-isotropic nature of the small- 
scale field. 

Upon further reflection on this problem, it is clear 
from the Biot-Savart law, which relates the veloc- 
ity field to the vorticity field, that anisotropy in the 
Reynolds stress can only arise if there is anisotropy 
in the two-point vorticity correlation tensor. Since 
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the anisotropy in Reynolds stresses persists at high 
Reynolds numbers, the vorticity field must also ex- 
hibit important anisotropy at high Reynolds numbers. 
The simulations led us to focus on this issue, and 
they suggest that modeling of the large-scale turbu- 
lence might be improved by explicit consideration of 
the anisotropies of the vorticity. 

Using the simulated fields as “data”, we explored a 
number of possible models for the return to isotropy 
problem, seeking a model that would describe the evo- 
lution more accurately. We found that the return- 
to-isotropy trajectories on the vorticity and dissipa- 
tion AIM’s were unique, and hence the evolution equa- 
tions for Vij and dij can be closed by simple models 
in terms of their own tensors. For example, the vor- 
ticity anisotropy was found to obey (during return to 
isotropy) 


dvjj 

dt 


-7-o v ij 

r 


( 10 ) 


Lee was able to obtain a tidy formula describing the 
coefficient 7 as a function of the invariants II b and II Vi 
and had similar success with the dissipation anisotropy 
model. Figures 3 and 4 show the comparisons of 
his model trajectories on the vorticity and dissipation 
AIMs for the re turn-to- isotropy cases. But unfortu- 
nately a truly satisfactory model for the evolution of 
b{j was not discovered, and so a complete second-order 
turbulence model has not yet been established. 

In summary, Lee’s simulations of homogeneous tur- 
bulence under irrotational strain provide some very 
important new insight into the physics of turbulence. 
They raise serious questions about traditional ideas of 
small-scale anisotropy that have been deeply imbed- 
ded in turbulence theory and modeling, and suggest 
that consider tat ion of small-scale anisotropy may be 
essential for development of better turbulence mod- 
els for the large-scale turbulence. Lee’s velocity fields 
have been archived at Ames and are being used to 
study other aspects of turbulence physics and model- 
ing. We expect that much more can be learned from 
them about turbulence. 


3. Scalar transport in homogeneous shear flow 

Homogeneous turbulent shear flow occurs when the 
mean velocity gradient transverse to the flow is uni- 
form. This situation has been studied experimentally 
in special wind tunnels that use variable grids to es- 
tablish this unusual mean flow, and has been explored 
in great depth using the Rogallo code. 

In homogeneous turbulence there is no net gain or 
loss to an elemental control volume as a result of tur- 
bulent transport (the gradients of all turbulence cor- 
relations vanish in homogeneous turbulence). But the 
turbulent fluxes are still present, and they can be stud- 
ied by direct simulation. Experiments on heat transfer 
in such a flow have been limited to cases with a mean 
temperature gradient aligned with the mean velocity 


gradient and cross-stream to it. These experiments 
suggest that there are strong anisotropies in the tur- 
bulent diffusion of scalars (e.g. temperature or species) 
in shear flow, but not enough is known about these ef- 
fects to allow them to be incorporated in turbulence 
models. 

In an effort to shed more light on this problem, Mike 
Rogers (1986), now a Staff Scientist at Ames, carried 
out detailed simulations of homogeneous shear flow 
with various arrangements for scalar transport, using 
the Rogallo code. The flow considered had the mean 
shear rate 

dlh 

dx 2 

and was examined with the three linearly-independent 
cases, each having a mean temperature gradient in 
only one direction, at Prandtl numbers of 1 , 0 . 2 , and 
0.7, and covered a wide range of Sq 2 /e, the dimen- 
sionless shear rate parameter, and Reynolds numbers. 
Most were carried out with a 128 3 mesh, in some cases 
non-square. The computed velocity and scalar fields 
have been preserved at selected time intervals, and rep- 
resent a great resource for subsequent study of new 
physical ideas or modeling concepts. 

The case of primary interest is case 2, having a mean 
temperature gradient in the x 2 direction, for this is the 
usual alignment of the velocity gradient and tempera- 
ture gradient in shear flows. The simulations confirm 
the rather surprising results of experiments that the 
heat flux in the direction perpendicular to the mean 
temperature gradient is significantly larger than the 
heat flux down the imposed mean temperature gra- 
dient! Hence, simple gradient-transport models that 
assume that the heat flux vector is aligned directly 
opposite to the mean temperature gradient will fail to 
predict the strong cross-gradient transport. There are 
many important practical situations where this coss- 
gradient transport is likely to be a very important phe- 
nomena, and so it is important that the mechanisms 
be understood and that the models be improved so 
that the fluxes can be predicted. 

Rogers and Moin studied the structure of the turbu- 
lence field in shear flow in great detail, using the simu- 
lations to explore such questions as the configurations 
of vortex lines in the flow. They found that hairpin 
vorticies, of the type found in inhomogeneous shear 
flows over solid walls, also occur in homogeneous shear 
flow. These vortices arise as a result of the straining 
of vortex filaments, and the legs of the hairpin vor- 
tices are predominantly aligned with the axis of pos- 
itive principal strain rate, i.e. at 45° upward from 
the x\ axis. In homogeneous shear flows there are 
equal numbers of “heads-up” and “heads-down” hair- 
pins. The vorticity in the heads has the same sense 
as the vorticity in the mean flow; indeed, much of the 
transverse vorticity is concentrated in these hairpin 
heads. Rogers also found that persistent shearing at 
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the higher Reynolds numbers gives rise to structures 
with strong spanwise coherence over substantial dis- 
tances, very much like the two-dimensional rollers in 
a two- stream mixing layer. Thus, the dominant large- 
scale structures in homogeneous shear flow are now 
thought to be some combinations of spanwise vortices 
and up and down hairpins. 

This model of the structure of the turbulence allows 
one to understand the anisotropic nature of the turbu- 
lent heat transport. Figure 5 shows Rogers’ explana- 
tion of the way in which these three types of coherent 
structures act to convect the scalar field in each of the 
three cases. 

Simple extensions of the gradient transport model 
can be used to predict the heat fluxes for arbitrary 
mean temperature gradients. In suggesting these mod- 
els we are not proposing that the transport process is 
one of gradient diffusion at the molecular level, but 
only that the mean temperature gradient provides the 
scaling of the heat flux. The direction of the heat flux 
vector is determined in a complicated way by the fluid 
dynamics, but it turns out that it is possible to model 
this in either of two rather simple ways. 

The first model uses an anisotropic tensor turbulent 
diffusivity Dij, defined by 

= (12) 

where 9* is the scalar(temperature) fluctuation, T is 
the mean scalar quantity (temperature), and F{ is the 
scalar flux in the i th direction. is easily calcu- 

lated from the simulation data, and is found to be 
non-symmetric. However, in a rotated coordinate sys- 
tem it can be made antisymmetric, and the important 
thing discovered by Rogers is that this coordinate sys- 
tem is essentially the same as the principal coordinate 
system of Rij and bij. This means that the diagonal 
members of Dij should be modeled in terms of the 
Reynolds stress tenor. In the rotated coordinate sys- 
tem, the off-diagonal elements must be antisymmet- 
ric, and the only antisymmetric tensor available is the 
mean rotation tensor, 


% = ( 13 ) 

This leads one to the model 

4 

D{j = Cibij -f C^ij + C^rUij) (14) 

where r is an appropriate turbulence time scale (q 2 /e) 
and the coefficients C\—Cz are dimensionless functions 
of scalars in the model. Rogers fit the coefficients as 
functions of the Reynolds and Prandtl numbers, and 
found that the resulting model could predict all of the 
heat fluxes for all of his fields to within 20%. This 


model can be used for inhomogeneous flows to estimate 
the anisotropic diffusivity tensor. 

A second, more sophisticated model, containing only 
one free coefficient, was also evaluated by Rogers. This 
model assumes that certain vector terms needing to 
be modeled in the exact equation for the turbulent 
scalar flux are aligned with the flux vector. This as- 
sumption had been made intuitively by earlier mod- 
lers, but Rogers was able to check it in detail using 
the simulations, and found that it was indeed approx- 
imately true. This assumption allows the equation for 
the scalar flux to be reduced to 





= 0 . 


(15) 


This is a set of linear algebraic equations for the Fi 
that can be inverted immediately for special gradients 
and solved locally in the general cases. Rogers fit the 
simulation results to determine C, and found that the 
resulting fluxes were predicted for all cases to within 
about 20%. 

In summary, Rogers’ simulations of scalar transport 
in homogeneous turbulent shear flow provide physi- 
cal insight on the structure of homogeneous turbulent 
shear flow, an explanation for the mechanisms of cross- 
gradient turbulent heat transfer, and simple models 
that can be incorporated into codes for more com- 
plicated engineering problems. The fields have been 
archived at Ames and continue to be explored in the 
study of the physics and modeling of turbulence. 

4. The spatially-developing mixing layer 

The mixing layer between two parallel streams of dif- 
ferent velocity has been widely studied in the labo- 
ratory and, for the case of time-development, by di- 
rect numerical simulation. The time-developing case 
has lent itself to simulation because it permits the use 
of periodic boundary conditions in the flow direction, 
thereby eliminating the need to devise turbulent inflow 
and outflow conditions. 

In a dissertation just completed, Pat Lowery (1986), 
who did his research as a NASA CFD Fellow and 
now works in CFD at the Battelle Northwest Re- 
search Laboratories, developed a code for the spatially- 
developing mixing layer, including scalar transport, 
with prescribed time-dependent inlet conditions. He 
studied the forced mixing layer in a series of two- and 
three-dimensional simulations, which included mixing 
studies and a fast chemical reaction calculation. 

Typical two-dimensional calculations used a mesh 256 
(cross-stream) by 1024 (streamwise). Three-dimen- 
sional calculations used a mesh 512 (streamwise) by 
128 (cross-stream) by 64 (spanwise). Finite-difference 
representations were used streamwise, spectral repre- 
sentations were used cross-stream and spanwise, and 
a third-order Runga-Kutta scheme was used for time 
advance. TVD corrections were necesary to maintain 
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proper scalar fields. On the Cray XMP 4-8, approxi- 
mately 60 CPU seconds were required per time step in 
three-dimensions. The Reynolds number at the inlet 
to the computational domain, based on the velocity 
difference across the mixing layer and the inlet vortic- 
ity thickness, was 100; the Prandtl number was 1. 

The first problem to be solved was the outflow bound- 
ary condition. Other workers have set gradients to 
zero at the outflow boundary, but this produces a 
“hard” exit condition that distorts the flow near the 
exit. Lowery found that a “softer” condition worked 
extremely well. The idea is simply to convect quan- 
tities out of the computational domain at a uniform 
convection velocity U c by setting 


d d 

* +v -arr° 


for all of the velocity componts and the scalar field. 
This allows the large-scale structures to roll out of 
the computational domain at the convection velocity, 
which was taken as the average of the two free-stream 
speeds. Figure 6 shows a snapshot of Lowery’s simu- 
lation, which strongly resembles pictures from experi- 
ments. 

Most of the work was done with forcing at the inlet, 
which in Lowery’s simulation was imagined to be just 
downstream of the splitter plate forming the mixing 
layer. Because the dynamical processes in the layer 
involve the interactions of modes with their subhar- 
monics to produce vortex merging (the mechanism by 
which the turbulence scale grows larger downstream), 
it was necessary to excite the inlet flow with fundamen- 
tals and sub-harmonics of the basic instability modes 
for the inlet laminar shear layer [tank profile). The re- 
sponse depends on the relative phases of these modes, 
and time permitted only a modest study of the pos- 
sibilities. However, it became clear (as experiments 
have already shown) that control over the growth of 
the mixing layer can be achieved by careful selection 
of the subharmonic phases. 

One of the major experimental observations is that 
more high-speed fluid is entrained into the layer than 
low-speed fluid. This effect is not captured in time- 
developing mixing layer simulations, since they are 
symmetric with respect to the two flow streams. This 
effect was captured in Lowery’s simulations, at a level 
commensurate with experiments. Probability density 
functions, the turbulent stress behavior, and other fea- 
tures of the flows were comparable in trends and mag- 
nitudes with experimental forced mixing layers, al- 
though it was impossible to make a direct comparison 
because of the sensitivity to the particular forcing. 

There are some theoretical arguments that suggest 
that the mixing layer is convectively unstable , mean- 
ing that removal of the forcing will permit the dis- 
turbances to wash downstream and the layer to re- 
laminarize. Lowery studied this question by freezing 


the disturbance field after some computation. Indeed, 
the disturbances washed downstream, but only part- 
way. It seems that the large-scale disturbances in the 
far-downstream region trigger instabilities in the up- 
stream flow through the pressure field. These new dis- 
turbances then travel downstream, repeating the trig- 
gering process (an absolute instability). It is generally 
believed that the splitter plate, absent in Lowery’s sim- 
ulation, is an essential element in making this flow ab- 
solutely unstable. The simulation suggests that this 
might not be the case. 

In his three-dimensional simulations, Lowery added 
streamwise vortices to the forcing. These produced 
structures very similar to those seen in the “braid” 
regions of the mixing layer. He ran one simulation 
with random excitation at the inlet, and found that 
these same structures were beginning to emerge natu- 
rally, although a great deal more computational time 
would have to be expended to allow them to develop 
completely. Thus, one can probably study these three- 
dimensional effects by judicious selection of the inlet 
forcing. 

Experiments have shown evidence of some vorticity 
of opposite sign to that in the mixing layer. This could 
be due to the vorticity coming from the boundary layer 
on the low-speed side of the splitter plate, or it could 
be due to overturning of vortex filaments by the three- 
dimensional action of the turbulence. Lowery could 
do what an experimentalist can not, i.e. run without 
the low-speed boundary layer, in which case negative 
vorticity could result only from overturning of vortex 
filaments. Feeding the inlet flow with vorticity of only 
one sign, only vorticity of this sign was found in the 
two-dimensional simulations. But in the three dimen- 
sional simulations a small amount of vorticity of oppo- 
site sign was found, indicating that the vortex filament 
overturning process indeed takes place. 

A new Ph.D. student, Neil Sandham, has used Low- 
ery’s code to study the mixing layer using inlet flows 
representing the splitter-plate boundary layers. These 
have shown very little effect of the boundary layers on 
the overall flow development, except very near the end 
of the splitter plate. In addition, we have developed a 
method for randomly exciting the inlet flow that looks 
quite promising in that it produces a mixing layer be- 
havior much more characteristic of natural mixing lay- 
ers than of the forced layers studied by Lowery. 

A new effort is being devoted to simulation of the 
compressible spatially-developing mixing layer, which 
is of interest in both aerodynamic and scramjet engine 
situations. The effort involves both direct numerical 
simulations at low Reynolds numbers and the develop- 
ment of improved models for LES studies of supersonic 
mixing and combustion. Supersonic mixing layers are 
known to spread and mix much less rapidly than sub- 
sonic layers, and the reasons for this are not well un- 
derstood. In the scramjet application, new ways to 
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enhance the mixing must be found. We are hopeful 
that these direct simulations will, over the next few 
years, shed some light on these important questions. 

In summary, Lowery’s simulations of the spatially- 
developing mixing layer provide new physical insight 
into the mechanisms of transport and entrainment in 
inhomogeneous free shear flows, and provide a basis 
for further studies on the modeling of such flows. In- 
let conditions for turbulent flows need to be developed, 
but the prospects look bright. Extensions to compress- 
ible flows are within reach and should be achieved in 
the next two or three years. 
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FIGURE l 

|Q V | v*. |n 6 I for relaxation from simple strains. • » AxR; A , BrR; 
V , CxR; O # DxR; □ , ExR; a # FxR; ■ f GlR; © , HxR where x = 
2, 3 , 4 . 4* , relaxation from axisymmetric contraction; x , relaxation 

from axisymmetric expansion. , locus of locking for relaxation 

from axisymmetric contraction (A$ = —1); — - , locus of locking for 
relaxation from axisymmetric expansion (A$ = +1). 
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Figure z 

lUrfl m. 1 14| for relaxation from simple strains. • , AxR; A , BxR; 
V , CxR; O , DxR; □ , ExR; a , FxR; a , GlR; O , HxR where r = 
2, 3, 4. +> , relaxation from axisymmetric contraction; x , relaxation 

from axisymmetric expansion. , locus of locking for relaxation 

from axisymmetric contraction (A* = —1); , locus of locking for 

relaxation from axisymmetric expansion (A* = +1). 


207 




C3R; A , C4R; V , E4R; O , F2R; □ , F3R; 9 , L6R; • , Q6R; ■ , 
L3W1R; , model prediction. 
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